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Abstract. We analyze in details the properties of the conditional state recently 
obtained by J. Wenger et al. [Phys. Rev. Lett. 92, 153601 (2004)] by means of 
inconclusive photon subtraction (IPS) on a squeezed vacuum state S(r)|0). The 
IPS process can be characterized by two parameters: the IPS transmissivity r 
and the photodetection quantum efficiency r\. We found that the conditional state 
approaches the squeezed Fock state >S(r)|l) when r, 77 — > 1, i.e., in the limit of singlc- 
00 1 photon subtraction. For non-unit IPS transmissivity and efficiency, the conditioned 

state remains close to the target state, i.e. shows a high fidelity for a wide range 
of experimental parameters. The nonclassicality of the conditional state is also 
investigated and a nonclassicality threshold on the IPS parameters is derived. 
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1. Introduction 

-i— > 

Beam splitters (BS) and avalanche photodetectors (APDs) play a fundamental role in 
Q-f quantum information processing. These key elements, among the other applications, 
can be used in order to generate non-Gaussian states from Gaussian ones [TJ El El and 
to distill continuous- variable entanglement [Sj. 

In this paper we focus our attention on the output state recently obtained 
experimentally by J. Wenger et al. [E] by means of photon subtraction on a squeezed 
vacuum state S(r) |0), S(r) being the squeezing operator. More precisely, when a 
Gaussian state, such as S(r) |0), is mixed with the vacuum at a beam splitter and, 
then, on/off photodetection is performed on the reflected beam, an unknown number of 
photons is subtracted from the input state and the output state is no longer Gaussian, 
i.e., the input state is de-Gaussified: this is due to the fact that the positive operator 
valued measure (POVM) describing the APD is non-Gaussian. Since the actual number 
of detected photons cannot be resolved by the APD, in [3] we referred to this process as 
to inconclusive photon subtraction (IPS). In general the IPS process can be characterized 
by two parameters: the beam splitter transmissivity r and the quantum efficiency 77 of 
the APD. As we will see, the conditional output state obtained by IPS on a squeezed 
vacuum is close to the squeezed Fock state S(r)\l), which is otherwise difficult to produce 
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Figure 1. Scheme of the IPS process: the input state g^ is mixed with the 
vacuum state go = |0)(0| at a beam splitter (BS) with transmissivity r; then, avalanche 
photodetection (APD) with quantum efficiency 77 is performed on the reflected beam. 
When the detector clicks we obtain the IPS state g( ont \ 



by Hamiltonian processes. For this reason, we address IPS as an effective resource to 
generate those squeezed Fock states. We find that the IPS conditional state reduces to 
S(r)\l) in the limit t, 77 — > 1, whereas for different values of the transmissivity and of the 
quantum efficiency it remains close to this target state, showing a high fidelity for a wide 
range of the parameters. Finally, since the IPS state obtained from the squeezed vacuum 
is, in general, non classical and mixed, we study how the purity and the nonclassical 
depth of the IPS state depend on r, 77, and on the input squeezing parameter r. 

The paper is structured as follows: in Section El we review the main elements of 
the IPS process on a single mode of radiation. The fidelity between the IPS conditional 
state and the squeezed Fock state S(z)\l), as well as its purity are then investigated in 
Section El whereas section |U is devoted to the analysis of the nonclassicality of the IPS 
state. Finally, Section closes the paper with some concluding remarks. 



2. The inconclusive photon subtraction process 

The scheme of the inconclusive photon subtraction (IPS) process is sketched in figure [D 
An input state is mixed with the vacuum state £0 — |0)(0| at a beam splitter (BS) 
with transmissivity r and, then, on/off avalanche photodetection (APD) with quantum 
efficiency 77 is performed on the reflected beam. Since the APD con only distinguish the 
presence from the absence of light, this measurement is inconclusive, namely does not 
resolve the number of the detected photons. In this way, when the detector clicks, an 
unknown number of photon is subtracted from the initial state and we obtain the IPS 
state g( ont \ Since the whole process is characterized by r and 77, we will refer to them 
also as IPS transmissivity and IPS quantum efficiency. 

If the input state of the mode a is the squeezed vacuum state Qr = |0, r) (0, r [, 
where |0,r) = 5(r)|0), S(r) = exp{^r(a^ 2 — a 2 )} being the squeezing operator (for the 
sake of the simplicity, without lack of generality, we can assume r as real), its (Gaussian) 
characteristic function Xr^(A a ) = x[^r m ^](A a ) reads 

xl in) (Aa)=exp{-iA>,.A a } (1) 
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where A = (x a , y a ) T , (• • -) T being the transposition operation, and 

cosh r + sinh r 

cosh r — sinh r 



(2) 



is the covariance matrix. Analogously, the vacuum state go = |0)(0| of the mode b is 
described by the (Gaussian) characteristic function 

Xo(Afc) = xM(A&) = exp{-|A^cr A b } , (3) 

where er = |l2, I2 being the 2x2 identity matrix. Since the initial two- mode state 
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Qo is Gaussian, under the action of the BS its 4 x 4 covariance matrix 
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transforms as follows [7j 
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is the symplectic transformation associated to the evolution operator of the BS. Now, 
the on/off photodetector with quantum efficiency n can be described by the POVM 
{n off (?7),n on (?7)}, with 



where A, B, and C are 2x2 matrices and 
Sbs = 



Vrl 2 


V1-t1 2 




v^l 2 



n fr(^) = $^(1 - r)) k \k)(k\, n on (77) = 1 - u oS ( v ) 



fc=0 



which corresponds to the characteristic functions 

x [n off (77)](A) =4 off) (A) = ~ex P {-§A T <7 M A} , 

X [n on (v)](A) = X JH(A) = 27n^(A) - X ( off )(A) , 
respectively, 5^ (A) being the 2-dim Dirac's delta function, and 
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The probability of a click in the detector is then given by [7] 
p on (r, r, 77) = Tr ab [g' r T I ® II oll (ry)] 

= / A a rf 2 A 6 X [< r ] (A , A 6 ) X [l] (- A a ) xJ on) (-A 6 

= 1- (*7 



'R 4 

v/DetfB + o-M] 



1 + (1 — Tgg) sinh 2 r 



-1 



(7) 

(8) 
(9) 

(10) 

(11) 
(12) 

(13) 



where x[^ T ](A a , A&) is the two-mode characteristic function associated to the state 
Q'rr = U BS £?! in) ® Qo C4 S , X[I](A) = 2tt<5( 2 )(A), and r cfr = r cff (r, 77) = 1 - 77(1 - r). Note 
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that when r c g — > 1, the probability (|13j) can be approximated at the first order in r c g 
as follows 

p on (r, t, rj) = (1 - r eff ) sinh 2 r + o[(l- r cfr ) 2 ] . (14) 
Finally, the output state 

conditioned to a click of the on/off photodetector, has the following characteristic 
function X r°r)j(A a ) = x[^r?](A ): 



(16) 



xS?(A») = o 7 \ / d 2 A bX [g' r , T ](A a ,A b ) X ^\-A b ) 

/ .[-r; °p{-Ks.M| (17) 

Pon(r,T,I,)\ 2 ' 1 «J + <T M ] J 

with Si = A and E 2 = A — C(S + ctm) _1 C t . Note that the output state is no longer a 
Gaussian state, namely its characteristic function is no longer Gaussian: for this reason 
the IPS process is also referred to as de-Gaussification process [§]. 

In general, a Gaussian state described by the characteristic function [in Cartesian 
notation, namely A = (x, y) T ] 

X (A)=exp{-±A T <rA} (18) 

with covariance matrix 



a c 
c b 



(19) 



can be also written in the complex notation as follows: 

X (\) = exp{-A\\\ 2 -B\ 2 -B*\* 2 } , (20) 

with 

A=±(a + b), B = i(b-a + 2k), (21) 

where we introduced the complex number A = 4=(x+iy). In this way, the characteristic 
function (|17|) can be written as follows: 

= expl-AIAI 2 -^ 2 -^* 2 } _ exp{-A 2 |A| 2 - ff 2 A 2 -g 2 *A* 2 } 

r,T,V Pon(r,r, V ) p on (r,T,r))r)^/Det[B + <T U ] ' 

where Ak and Bk are refers to the covariance matrix S^, k = 1,2 respectively. Finally, 
using the definition 

W[g](a) = ^- / d 2 \x[g](\) exp{A*a-a*A} , (23) 
Jc 

which relates the Wigner function W[g](a) of a state g to its characteristic function 
x[^](A), one can obtain the Wigner function Wr°^(a) = W[^r°^](cK). As for the 
characteristic function, to pass from the complex, W[g](a), to the Cartesian notation, 
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Figure 2. (a) Plot of the Wigner function W}°t, v (x,y) with r = 0.5, r = 0.90, 
and r\ — 0.80; (b) plot of the Wigner function Wz SqF \x,y) of the state S(z)\l) with 
squeezing parameter z = 0.5. 



W[g](x,y), one should put a = -^{x + iy) [Zj- In figure El (a) we report Wr°^ (x,y) for 
fixed r, r, and n: as it is apparent from the plot the Wigner function is not Gaussian, 
and may assume negative values jO] . In Section |U we will investigate this effect by 
analyzing the nonclassicality of the conditioned state. In figureEl(b) we show the Wigner 
function x* ( x > v) associated to the squeezed Fock state = S(z)\l)(l\S^(z), 

whose characteristic function x^ (A) = x[Qz ](A) reads (we assume z as real) 

xi SqF) (A) = [1 - 2 (A|A| 2 + £ A 2 + B* \* 2 )} exp {-A|A| 2 - B X 2 - B* X* 2 } , (24) 

with Ao = 2(cosh 2 2; + sinh 2 2;) and Bo = —2 cosh z sinh^;. Since the Wigner functions 
of the IPS squeezed vacuum and of the squeezed number state are quite similar, one 
can think of using the IPS process to produce the state ^r SqF ' ) ; motivated by this 
consideration, in the next section we will analyses the fidelity between this states. Figure 
E] shows Wr°rtf {x,y) with fixed r and rj and different values of the IPS transmissivity 
r; the plots on the right of the same figure compare the W^rn (0, y) (solid lines) with 
Wr SqF \o, y) (dashed line). Finally, the effect of the quantum efficiency rj on the output 
state is shown in figure 0] where we plot as reference the value of the Wigner function 
Wr.T^ at the center of the complex plane as a function of the transmissivity r and 
different values of n: we can see that the main effect on the output state is due to r. 
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Figure 3. Plots of the Wigner function Wr™ v (x,y) with r = 0.5, -q = 0.80 and 
different values of the BS transmissivity r: from top to bottom r = 0.99, 0.9, 0.75, 
and 0.50. The solid lines of the plots on the right refer to Wr°r^ (0, y) whereas the 
dashed lines are Wz (0,y) of the state S(z)\l) with squeezing parameter z = 0.5. 
Note that when r = 0.99 the two lines overlap, y is the squeezed coordinate. 



3. Fidelity and purity 

The fidelity between the pure state £?i SqF ^ and the IPS state £r°"r? is defined as follows: 

F T , v (z,r) = Tr[ g ^^] (25) 

= f / rf 2 Axi SqF) (A)x^S(- A )' ( 26 ) 
2tt J R 2 



where 



Pon(r,r,T]) 1 rjy/Det[B + cr M ] 



= Al-Al- 4(gg - Bp 
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Figure 4. Plots of Wr°T,r, / (0, 0) with (a) r = 0.5 and (b) r = 2.0 as a function of r 
and different values of r\: from bottom to top r\ — 1.0, 0.75, 0.50, and 0.25. The value 
of the function is mainly affected by t. 
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Figure 5. Plot of the fidelity F TtV (r) with rj = 0.80 as a function of the IPS 
transmissivity r for different values of r: from top to bottom r = 0.1, 0.3, 0.5, 0.7, 1.0, 
and 2.0. 



and Ah and B^, /i = 0, 1, 2, have been introduced in equations (jZH) and (|2*2*j) . respectively. 
The analytic expression of F TjV (z,r) is quite cumbersome, but, on the other hand, we 
can draw some interesting consideration by addressing its expansion at the first order 



in the transmissivity r when r 
F T Az,r) ' 



1 and rj — 1, namely 



cosh (r — z) 
9 cosh(r + z) 



3 cosh(3r 



;i-r) + o[(l-r) 2 ] . 



(29) 



8 cosh(r — z) 

In fact, from the expansion ()29|) we conclude that the maximum of the fidelity is achieved 
when z = r. 
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Figure 6. Plot of the purity fi T . v (r) of the state Q °T,-q- We set 77 = 0.80. 



In figure El we plot F TtV (r) = F TtV (r,r) as a function of the IPS transmissivity 
and for different values of r. We can see that F T<n reaches it maximum when the IPS 
transmissivity approaches 1, namely in the single-photon subtraction limit j3]. Moreover, 
when the squeezing parameter r increases the fidelity decreases: this is due to the 
increasing (unknown) number of subtracted photons which reduces the purity of the 
IPS state itself. In figured we plot the purity /i Ti r ? ( r ) of the IPS squeezed vacuum q^t^ , 
defined as follows [Bj: 

^(r) = Tr [{q%$)*\ = vr / d 2 a [W$${a)] 2 (30) 

1 r 1 1 



2p on (r,T,ri) \ y/Al-ABf v 2 Bet[B + a M } \J A\- AB\ 

4 

~ V ^Det[B + ^(A + A 2 ) 2 - Wi + W 



(31) 



4. Nonclassicality of the IPS squeezed vacuum state 

As a measure of nonclassicality of the IPS state gi^rj we consider the nonclassical depth 

m 

r=^, (32) 
s being the maximum s for which the generalized quasi-probability function 

W s (a) = - [ d 2 \ X (\) exp{is + X*a-a*\} (33) 

is a probability distribution, i.e. positive semidefinite and non singular. As a matter 
of fact, one has T = 1 for number states and T = for coherent states. Moreover, 
the nonclassical depth can be interpreted as the minimum number of thermal photons 
which has to be added to a quantum state in order to erase all the quantum features of 
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the state El] . In the case of gi^rj , we have [for the sake of simplicity we do not write 
explicitly the dependence on r,r and n in the symbol wj out ^(a)] 

ut) («) = / r UM) ~ , GM) ) (34) 

where we defined 

f 2(2A k - s)\a\ 2 + 4Bja 2 + 4B k a* 2 
g() 6XP l (2A k - s) 2 - 16\B k \ 2 

ny/(2A k -sy- 16\B k \ 2 

At first we note that in order to have Ws (a) normalizable, equation the following 
condition should be satisfied 

s < 2A k (k = 1,2). (36) 

Furthermore, since is a difference between two Gaussian functions with the 

center in the origin of the complex plane, one can easily see that, in general, this function 
has a minimum in a = and that this minimum can be negative. For this reason and 
thanks to other simple considerations about the symmetries of wj ont \a) with respect to 
the point a = 0, we can focus our attention in the origin of the complex plane, obtaining 
this further condition for the positivity: 

^(0)-^=§|L=>0, (37) 

which, together with the conditions ([36)1. brings to 

2 - 7] - (4 - V )t 

S ^ V) = 2-(l-r), ' (38) 
and, then, to the following expression for the nonclassical depth: 

T(t,v) = „ « x • (39) 
2-(1-t)t7 

Since T(r, rf) > 0, the conditional state is nonclassical for any non-zero value of the IPS 
transmissivity and efficiency Note that equation ()39|) depends only on r and r], whereas 
it is independent on the squeezing parameter r. Notice, however, the nonclassical depth 
does not measure the extension of the negativity region, but only the presence of negative 
values. Therefore it is not surprising that equation (|3U|) does not depend on r. We 
plot T(t, if) in figure [7| Since the usual Wigner function is obtained when s = in 
(jHSl) , from equation (J3~%|) we can see that W^S (a) becomes semi-positive definite when 
r= (2 - 77)/(4 - 77). 



5. Concluding remarks 

We have analyzed in details the state obtained subtracting photons from the squeezed 
vacuum by means of linear optics, namely using beam splitters and avalanche 
photodetectors. We referred to the whole photon-subtraction process as to inconclusive 
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Figure 7. Plot of the nonclassical depth T(r, rf) of the IPS squeezed vacuum state. 

photon subtraction (IPS), since avalanche photodetectors are not able to resolve the 
number of detected photons. We found that the IPS conditional state obtained from a 
squeezed vacuum state is close to the squeezed Fock state S(r)\l) and approaches this 
target state when only one photon is subtracted, namely, using a high transmissivity 
beam splitter for the IPS. Moreover, when the transmissivity and the quantum efficiency 
are not unitary, the output state remains close to the target state, showing a high 
fidelity for a wide range of the parameters. The purity and the nonclassicality of the 
IPS squeezed vacuum state have been also considered: we found that the relevant 
parameter is the transmissivity r, while the IPS efficiency n only slightly affects the 
output state. We conclude that IPS, which was recently experimentally implemented 
jH], can be effectively used to produce a nonclassical state such as the squeezed Fock 
state S(r)\l), whose generation would be, otherwise, quite challenging. 
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